AO-A044  912 


UNCLASSIFIED 


AO 

A0449I2 


STANFORD  UNIV  CALIF  DEPT  OF  OPERATIONS  RESEARCH  F/6  l2/i 

A LEAST  ELEMENT  THEORY  OF  SOLVINO  LINEAR  COMPLEMENTARITY  PR0BLE~ETC(U) 
dec  76  R M COTTLE » J PANS  N00014-75-C-0267 

TR-76-30  NL 


A LEAST  ELEMENT  THEORY 


OF  SOLVING  LINEAR  COMPLEMENTARITY  PROBLEMS 
AS  LINEAR  PROGRAMS 

by 

Richard  W.  Cottle  and  Jong-Shi  Pang 


TECHNICAL  REPORT  76-30 
December  1976 


DEPARTMENT  OF  OPERATIONS  RESEARCH 

Stanford  University 
Stanford,  California 

Research  and  reproduction  of  this  report  were  partially  supported 
by  National  Science  Foundation  Grant  MCS  71-03341  A04  and  Office 
of  Naval  Research  Contract  N00014-75-C-0267 . 


Reproduction  in  whole  or  in  part  is  permitted  for  any  purposes 
of  the  United  States  Government.  This  document  has  been  approved 
for  public  release  and  sale;  its  distribution  is  unlimited. 


A LEAST  ELEMENT  THEORY 
OF  SOLVING  LINEAR  COMPLEMENTARITY  PROBLEMS 
AS  LINEAR  PROGRAMS 

by 

Richard  W.  Cottle  and  Jong-Shi  Pang 
ABSTRACT 

In  a previous  report  [2],  the  authors  have  established  a 
least-element  interpretation  to  Mangasarian ' s theory  [5],  [6]  of 
formulating  some  linear  complementarity  problems  as  linear  programs. 
In  the  present  report,  we  extend  our  previous  analysis  to  a more 
general  class  of  linear  complementarity  pr  ’-''ems  investigated  in 
Mangasarian  [7].  Our  purposes  are  (1)  t onstrate  how  solutions 

to  these  problems  can  be  generated  from  least  elements  of  polyhedral 
sets  and  (2)  to  investigate  how  these  "least-element  solutions"  are 
related  to  the  solutions  obtained  by  the  linear  programming  approach 
as  proposed  by  Mangasarian.  1 


A LEAST  ELEMENT  THEvORY 


OF  SOLVING  LINEAR  COMPLEMENTARITY  PROBLEMS 
AS  LINEAR  PROGRAMS 

by 

Richard  W.  Cottle  and  Jong-Shi  Pang 
1 . INTRODUCTION 

In  this  paper,  we  study  the  linear  complementarity  problem  of 
finding  a vector  x e R*^  satisfying 

(1.1)  X ^ 0,  q + Mx  > 0 and  x"^(q  + Mx)  = 0 

where  the  given  n-vector  q and  n x n matrix  M satisfy  the 

following  three  assumptions: 

(Ml)  MX  = Y + qc’^ 

(M2)  r"^X  + s'^Y  > 0 

(M3)  r'^(X+C)  + s^(Y-K:)  > 0 

whre  X and  Y are  suitable  Z-matrices  (i.e.,  real  square  matrices 
whose  off-diagonal  entries  are  non-positive) , C is  a diagonal 
matrix  whose  diagonal  elements  are  the  components  of  the  vector  c, 
and  r,  s and  c are  some  non-negative  vectors.  We  denote  problem 
(1.1)  by  the  pair  (q,M).  Its  feasible  set  is  defined  as  the 
polyhedral  set 

X(q,M)  = {x  e R*’  : x > 0,  q+Mx  ^0' 


1 


The  linear  complementarity  problem  (q,  M;  with  the  vector  q 
and  matrix  M satisfying  conditions  (Ml)  - (M3)  has  recently  been 
studied  by  Mangasarian  [?]  who  shows  that  such  a problem  can  be 
formulated  as  the  linear  program 

T 

(1.2)  minimize  px  subject  to  x>0  and  q + Mx  > 0 

T 

where  p = r + M s.  Our  purpose  in  this  report  is  to  show  that  this  result 
is  related  to  a theory  of  polyhedral  sets  having  least  elements.  (A  vector 
X belonging  to  a set  S C is  said  to  be  the  least  (greatest)  element 
of  S if  X < ( >) -X  for  every  x t S.)  The  method  of  derivation  used 
by  Mangasarian  is  not  based  on  least-element  arguments.  In  a previous 
report  [2j,  the  authors  have  applied  this  theory  of  polyhedral  sets  having 
least  elem.ents  to  the  particular  case  c = 0 and  established  that  for 
every  n-vector  q,  the  linear  complementarity  problem  (q,  M)  has  a 
solution  which  can  be  generated  from  the  least  element  of  a polyhedral  set, 
thus  providing  a least-element  interpretation  to  the  linear  programming 
formulation  of  the  problem  (q,  M)  which  was  initially  obtained  by 
Mangasarian  in  [5].  In  the  present  report,  we  extend  our  previous  analysis 
to  the  general  case  where  c is  merely  non-negative,  as  described  at 
the  beginning  of  the  introduction.  Our  purposes  are  (l)  to  establish  the 
least-element  characterization  of  a solution  to  the  linear  complementarity 
problem  under  consideration  and  (2)  to  demonstrate  how  this  "least-element 
solution"  is  related  to  the  solution(s)  obtained  by  the  linear  programming 
approach  as  proposed  by  Mangasarian.  Here,  we  should  point  out  that  all 


■ 

the  lliicor  complc'incnUiiity  problem'.,  whtcli  ere  listed  In  Teble  1 of  [ 7 ) to  be 
solvoble  es  linear  progra/ns,  satisfy  conditions  (Ml)  - (MJ).  Therefore  the 
least-element  theory  developed  In  the  present  report  is  applicable  to  all  of 
them.  However,  it  may  not  necessarily  be  applicable  to  those  problems 
satisfying  the  more  general  conditions  in  Theorem  1 of  [?). 

It  would  bo  appropriate  for  us  to  review  some  of  the  essential  results 
obtained  in  [2]  for  the  particular  case  c = 0.  Using  the  same  notations, 
we  denote  problem  (1.2)  by  the  triple  (p,Q,M)  and  by  C the  class  of 
square  matrices  M for  which  there  exist  Z-matrices  X and  Y such  that 
the  following  two  conditions  are  satisfied 


(Ci)  MX  = Y 

T T 

(C2)  r X + s Y > 0 for  some  r,  s > 0 . 

These  are  precisely  conditions  (Ml)  - (M3)  v/ith  c = 0.  The  following 
proposition  is  an  immediate  consequence  of  the  well-known  theorem  of 
Kuhn-Fourier  [4]  on  the  solvability  of  a system  of  linear  relations. 
Proposition  1 . I . Lot  X and  Y be  n Xn  matrices.  Then  the  following 
are  equivalent 


(C2) 


(C2)' 


T 

r X 


T 

s Y > 0 


for  some  r,  s > 0 


u > 0 


Xu  < 0 > =>  u = 0 . 

Yu  < 0 
— 


We  have  established  useful  necessary  and  sufficient  conditions  for 
two  Z-rnatriccs  X and  Y to  satisfy  condition  (C2).  These  are  stated  below. 


I 

< 

i 

1 

I 

I 
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I.c‘miiui  )■?.  L«,‘t  X and  Y bo  X-mctriccs  of  the  fiomo  order.  Tlu;n  (C2) 

holds  if  and  only  if  there  exist  a principal  rearranyeinent  with  permutation 
matrix  P and  a partitioning  of  X and  Y such  that 

(1.3a)  P^XP 


(y  Y \ 

11  12 

T 

, P YP  = 

111  12 

X.,  X’  , J 

Y Y 

21  ZZj 

1 21  22/ 

( 1 • 3b) 


X X ^ 
11  12 


Y Y 
21  22 


is  a K-matrix 


Using  this  lemma,  wo  have  given  necessary  and  sufficient  conditions 
for  M t C . 

Theorem  1-3.  Let  M,  X and  Y be  n X n matrices  with  X and  Y 
both  2-m.atrices.  Then 
(Cl)  MX  = Y 

(C2) 


T T 

r X + s Y > 0 for  seme  r,  s > 0 


if  and  only  if  there  is  a principal  rearrangement  and  partitioning  of  M,  X 
and  Y such  that 


“l2\ 

/^1 

^2l 

^2' 

(l.4d) 

i“21 

“22/ 

t^21 

\z\ 

i"2. 

''22; 

f^i 

^2l 

(1.4b) 

i’'2, 

1 is  a 

K-matrix 

(1.4c) 

X 

is 

nonsingular  . 

A K-matrix.  also  known  as  Minkc>wski  matrix.  Is  a Z-i:iatrix  with  a non- 
negative Itivorsc. 
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The  above-mcMiUoned  leiist-eleincnl  resell  foi  Ihe  linear  compleit.enlui ity 


problem  (q,  M)  with  M t C is  stated  in  the  theorem  bolcw. 

Theorem  1.1.  Let  M t C ond  let  X and  Y bo  Z-malrices  satisfying 
(Cl)  and  (C2).  Suppose  (q,  M)  is  feasible,  i.e.  X(q,  M)  / 4>.  Then 
the  polyhedral  set 

V = {vt  R*^:Xv^0,  q+Yv>0} 

contains  a least  element  v.  Moreover,  the  vector  x = Xv  solves  the 
problem  (q,  M). 

As  a consequence  to  this  theorem,  we  deduced  that  for  every  vector 
q t the  linear  complementarity  problem  (q,  M)  with  M c C can  be 

solved  as  the  linear  program  (p,  q,  M)  where  the  vector  p is  the  (unique) 
solution  to  the  system  of  equations 

T 

p X = f 

for  some  positive  vector  f.  We  have  also  shown  that  the  vector  p 
required  in  Mangasarian's  theory  can  be  obtained  in  precisely  the  same 
way.  In  the  last  part  of  the  report,  we  established  several  related  matrix- 
theoretic  results,  and  demonstrated  that  C includes  all  the  matrices 
investigated  by  Mangasarian  in  [6]. 

We  explain  the  notations  used  in  the  paper.  All  vectors  and  matrices 
under  consideration  are  real.  A Z-matrix  X is  said  to  be  a ^-niQlnx 
if  (X  + el)  is  a K-matrix  for  every  e > 0.  The  letters  Z,  and  K 
will  also  denote  the  class  of  Z-,  K^-  and  K-matrices  respectively. 
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Various  characterizations  of  K-  and  KQ-inati ices  can  be  found  in  I 3). 
Wo  denote  the  range  space  of  a matrix  A by  W(A),  i.e.  ft(A)  consists 
of  those  vectors  which  can  be  represented  as  linear  combinations  of  the 
columns  of  A.  Let  M be  an  n x n matrix.  If  I,  J L • • • > n}, 
we  define 


m.  . 

Vi 


m.  , 
‘s^ 


m. 


m.  , 
^s^ 


where  1=  {i,...,i  } and  J = {j  i.)  with  l<i,  <•••<  i 

IS  It  la 

and  1 < j j < • • • < jj.  < n.  In  particular,  is  a principal  submatrix 

T 

of  M Similarly,  if  q c R , we  define  = (q^  , .,q.  ) . We  denote 


th 


ation  vector  (l,...,l)  by  e. 


2 . CONNlxrnON  WITH  l.LAST  lii.ll.MIiN'r.S 

2.1.  Gcnoigl  di.scussicm.  Lot  q bo  on  n • vector  and  M an  n X n 

matri.x  satisfying  conditions  (Ml)  - (M3)  for  some  Z-matriecs  X and  Y, 

and  some  non-negativo  vectors  r,  s and  c.  Throughc;Ut  this  section, 

these  vectors  and  matrices  are  assumed  to  possess  the  properties  just 

mentioned.  We  shall  develop  a least-element  study  of  the  linear 

complementarity  problem  (q,  M)  with  such  a vector  q and  matrix  M. 

As  a consequence  of  our  investigation,  we  shall  establish  a least-element 

interpretation  for  the  result  obtained  by  Mangasarian  [7  ] of  formulating 

such  a linear  complem.entarity  problem  as  the  linear  program  (p,  q,  M)  with 
T 

p = r + M s.  We  start  by  proving  a lemma  which  strengthens  condition  (M3). 
Lemma  2.1.  Let  c,  r and  s be  non-negative  vectors  and  let  X and 
Y be  Z-matrices.  If  condition  (M3)  holds,  then 
(M3)'  r + s > 0 . 

Proof:  It  suffices  to  show  that  for  every  i = 1,  . . . , n,  r.  = 0 (s.  = 0)  =>  s.  > 0 
(r.  > 0).  So  assume  r^  = 0,  say.  Then 

0 < (r^(X  + C)  + s'^(Y  + O). 

= Tj  r.X..  + X s.Y  . + s.(Y.,  + c.) 

. j.  1 Ji  “T-  1 11  1 11  1 

jyi  )i-i  ' ^ 

< s ,(Y, . + c.)  . 

~ 1 11  1 


Thus  Sj  > 0.  Similarly,  we  may  deduce  [ = 0 =>  r.  > 0].  Therefore 

(M3)'  follows.  n 
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1 

R(Miuirt. . CoiKiillon  (Mi)'  was  rcluirod  to  os  a spriciiil  case,  biii  not  j.*;  o 
consocjuodcc  oi  (M  i)  in  1 V ] . In  fact,  if  c > 0 (as  in  Cor(<H.iry  1 of  ( 7 |), 
the  two  conditions  arc  equivalent. 

Wo  recall  that  if  M t C,  the  matrix  X salisfyinn  (Cl)  and  (C2)  j 

j 

must  bo  nonsingular  (Theorem  1.3).  The  following  example  illustrates  ^ 

that  there  can  be  singular  matrix  X satisfying  (Ml)  - (M3). 

Exam  [lie  P..Z.  Let  q = (_^)  and  M = (J  _J).  The  problem  'q,  M;  has 
two  solutions,  namely,  (j/p  ^ ^ 0 0^’  ^ ~ ^0  -p’ 

r = (^)^  s = (j',)  and  c = (1),  then  conditions  (Ml)  - (M3)  are  satisfied. 

Nevertheless,  X is  singular.  M.ireover,  it  is  not  hard  to  verify  that  M {C.  . ] i 

The  fact  that  there  exist  such  singular  matrices  X indicates  that 
in  order  to  develop  a least-element  theory  for  the  linear  complementarity 
ifrofiTenf^^q,'^/) , one  should  not  merely  concentrate  on  the  range  space 
of  X.  In  fact,  the  same  example  above  shows  that  R(X)  fl  X(q,  M)  = <(). 

Later  in  our  discussion,  we  will  see  that  this  latter  relation  always  holds 
if  X is  singular  and  the  stronger  condition  (C2)  is  imposed  (see 
Proposition  2.26  and  the  remark  following  it). 

It  is  clear  that  if  M c C,  then  conditions  (Ml)  - (M3)  are  satisfied 
for  every  vector  q.  Nevertheless,  if  a matrix  M satisfies  (Ml)  - (M3) 
for  some  vector  q,  it  does  not  necessarily  follow  that  M c C.  Example  2.2 
illustrates  this  fact.  The  following  provides  another  example. 
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1 1, 


Excuiiijlc  3,.. 3.  Lot  q = (_^)  and  M = ( ^ _j).  Then  conditiimb  (Ml)  - (M3) 
are  satisfied  with  r = (^),  s - c = X = (~J  and  Y = ( f °). 


We  show  M ^ C.  Suppose  not,  then  there  exists 


' ^11  "^1? 
"^21  ^22 


1 0' 
t Z such  that 


'-1 

1\ 

X 

1 1 

-"12^ 

■’‘21 

X 

12 

-22  \ 

2 

V 

-M 

^"''21 

^zzi 

"‘21 

-2x  ^ 
12 

-J 

« Z 


and  (r 


•■v'  . 


1 1 

-^2 

21 

^22 

21 


^‘21  ^ ^22 


■n  -"hz-'zz 


> 0 for 


some 


r 2'  -X. 

' ^ 

non-negative  scalars  r ,r.,s,  and  s^.  It  follows  that  2x  + x,  < 0 

1’  2’  1 2 1121" 


which  implies  ^j2-"2^21-°'  i^^plies 

I-emma  1. 2 implies  that 


,2x^1 


e K . 


In  particular,  det 


-X  - x^ 
11  21 

2x  + X,, 
11  21 


^12  ^22 
■"’‘12  -’‘22/ 


> 0.  Since  det  M < 0,  it 


follows  that 


(2.  I) 


det 


^1 

"^21 

X 1 
"^22/ 

= ’‘u’‘22-’‘l2’‘21  ■ 


Furthermore, 


We  have  shown  that  < -x  < 0;  thus  >=,,Xp,>x  x contradicting 


(2.1).  Therefore  ( * ^ C. 


22 
-1  1, 
2 


n 
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The  obovc  examples  illustrate  that  the  present  class  of  linear 
complementarity  problems  (q,  M)  with  q and  M satlsfyiny  (Ml)  - (M3) 
is  a genuine  extension  of  the  previous  class  of  problems  (q,  M)  with 
M t C . Example  2.  3 was  used  or'ginaliy  in  [ 7 ] by  Mangasarian  for 
another  purpose. 

2.2.  The  set  "U . V/e  define 

V.  = {u  t ; u > 0,  Xu  < 0,  Yu  < 0,  c^u  = 1 } . 

This  polyhedral  set  l(  plays  a very  important  role  throughout  our  whole 
discussion.  It  may  be  empty,  for  example,  if  c = 0 as  in  the  case  of 
M € C . The  following  example  illustrates  that  it  can  sometimes  be  non- 
empty as  well. 


0 

3 

4.  ,-2 

1 

-1 

0 and  q = 0 . Then  v/ith  c 

0 

-1 

-3  1 \ 1 1 

2 

-1 

-2 

-1 

1 

0 , conditions  (Ml)  - (M3)  are 

-1 

0 

2 

<2/5 

satisfied.  In  this  case,  u = 2/5  € T<.  We  shall  say  more  about  this 

W/si 

example  later.  □ 

The  proposition  below  describes  the  relationship  between  vectors 
in  V (if  any)  and  solutions  to  the  linear  complementarity  problem  (q,M). 
Proposition  2.5.  If  u < then  x = -Xu  is  a solution  to  (q,  M). 

Proof:  If  X = -Xu  where  u r Tf,  then  x > 0;  moreover, 
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q + Mx  --  q - MXu 


= q - (Y  + qc')u 
= -Yu  > 0 , 

i.c.  xtX(q,M).  We  also  have 

T T 

0<rx  + s (q  + Mx) 

- -(r'^X  + s^Y)u  < 0 . 

Therefore,  for  each  i = 1,  . . • . , n, 

rfXf  = s.(q  + Mx),  = 0 . 

If  X,  > 0,  say,  then  r,  = 0.  Condition  (M3)'  implies  s,  > 0 which 
gives  (q  + Mx).  = 0.  Therefore  the  vector  x defined  above  solves  the 
problem  (q,  M).  □ 

Remark  I.  In  the  proof  above,  the  assumption  that  X and  Y are  Z-matrices 
is  required  in  order  for  condition  (M3)'  to  be  applicable.  Therefore,  if 
condition  (M3)'  holds  by  itself  (as  in  Corollary  1 of  [ 7 ]),  then  the 
proposition  is  valid  for  any  matrices  X and  Y (which  are  not  necessarily 
Z-matrices)  satisfying  conditions  (Ml)  and  (M2). 

Remark  2.  The  feasibility  of  the  problem  (q,  M)  is  not  a requirement, 
but  a consequence  of  Proposition  2.  5. 

Proposition  2.  5 shows  that  if  the  set  is  nonempty,  then  a 
solution  to  the  linear  complementarity  problem  (q,  M)  can  be  obtained 
by  first  finding  any  vector  u in  K , or  equivalently,  solving  the  system 

of  linear  inequalities 
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u > 0,  Xu  < 0,  Yu  < 0 und  c u - I 

and  thcMi  setting  x = -Xu.  The  question  of  whether  is  noncMiipty  is 
itself  answerable  by  linear  progiamniincj.  The  next  proposition  shows 


that  if  T/  ^ 4>,  then  any  vector  x such  that  x = -Xu  for  some  u € 


actually  solves  the  linear  program  (p,q,M)  where  p = r + M s. 


Proposition  2. 6.  If  u t then  x = -Xu  solves  the  linear  program 

.T 


(p,  q,M)  with  p = r +M  s. 
Proof:  If  X e X(q,  M),  then 


TXT  T 

p X = (r  4 s M)x  ^ -s  q . 


On  the  other  hand,  the  proof  of  the  last  proposition  shows  that 
X = -Xu  t X(q,  M)  if  u t "U  . Furthermore,  with  such  an  x,  we  have 


T T T 

p X = -(r‘  + s M)Xu" 


, , T„.-  T T 

= -(r  X 4-  s Y)u  - s q = -s  q . 


Therefore,  any  such  x solves  the  linear  program  (p,  q,.M).  □ 

Corollary  2.7.  If  M is  nondegenerate  (i.e.  every  principal  submatrix 

of  M is  nonsingular)  and  if  V * then  the  linear  program  (p,  q,  M) 

T — — 

with  p = r 4-  M s has  a unique  solution  x.  Furthermore,  x = -Xu  for 

every  u t %( . 

Proof;  That  the  linear  program  has  a solution  x follows  from  Proposition  2.6. 
To  show  its  uniqueness,  let  x be  a solution  to  (p,q,M).  Then  the 
proof  of  Proposition  2.  6 shows  that,  in  fact. 


r^x  = s^(q  4-  Mx)  = 0 
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If  I = {l  : > 0}  and  J = (1,  . • then  x must  necessarily 

satisfy  the  following  system  of  equations 

( 0 \ f 1 0 

(2.2)  + 

\qjl  \ Mj, 

Here  we  have  used  the  fact  that  s^  > 0 by  (M3)'.  By  the  nondegeneracy 
of  M,  the  system  (2.2)  ha.s  a unique  solution.  This  establishes  the 
uniqueness  of  x.  The  last  conclusion  of  the  corollary  is  immediate.  ,3 

Remark.  The  matrix  M in  Example  2,4  is  degenerate:  nevertheless,  the 
proof  of  the  corollary  shows  that  any  solution  to  the  linear  orogram  (p,  q,  .M) 
must  satisfy  q + Mx  -.0.  Since  M is  nonsingular,  the  LP-soiution 


X 


X 


2/5, 

2/5 

1/5 


is  unique.  This  solution  was  also  obtained  in  [ 7 ] by  actually 


applying  the  simplex  method  to  the  linear  program.  Its  uniqueness  was 
not  proven  there. 

The  results  above  establish  the  relationship  between  vectors  in  V, 
solutions  to  the  linear  complementarity  problem  (q,  M)  and  solutions  to 
the  linear  program  (p,q,  M),  In  the  sequel,  we  divide  our  discussion  into 
cases.  In  each  case,  we  shall  construct  a polyhedral  set  having  a least 
(or  greatest)  element  and  demonstrate  how  a solution  to  the  linear  comple- 
mentarity problem  can  be  generated  from  this  element.  In  the  analysis,  we 
shall  need  two  fundamental  results  from  lattice  theory.  These  are  stated 
in  Propositions  2.9  and  2.10.  Their  proofs  are  easy  and  are  omitted. 
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Definition  2.  K.  Let  S be  n subset  of  Then  S Is  said  to  be  a 

meet  (join)  r.ornl-sublnttico  (of  H*^)  if  for  every  s,  t c S,  the  meet  ( join) 
of  s and  t,  defined  as  the  vector  u - (u^)  where  - min(s.,t.) 
(max(s, , t.))  for  every  i,  belongs  to  S.  The  set  S is  bounded  below 
(above)  if  there  exists  a vector  s e R such  that  s > (<)  s for  every 

s e S. 

Proposition  Z.  9 . The  following  are  equivalent: 

1)  I.  is  a polyhedral  meet  (join)  semi-sublattice  of  R^. 

2)  L = {s  « r"^  : As  > b}  for  some  matrix  A and  vector  b,  with  A 
having  at  most  one  positive  (negative)  element  in  each  row. 

Proposition  2.10.  Let  L be  a nonempty  meet  (join)  semi-sublattice  of  R*^. 
If  L is  closed  and  bounded  below  (above),  then  L has  a least  (greatest) 
element.  Further.more,  this  element  can  be  obtained  by  solving 

T 

(2.3)  minimize  (maximize)  f x subject  to  x « L 

for  any  positive  vector  f. 

We  should  also  mention  that  a least  (greatest)  element  of  a meet  (join) 
semi-sublattice  L of  R*^  actually  solves  the  program  (2.  3)  for  any  non- 
negative vector  f and  it  is  the  unique  solution  if  f is  positive. 

In  the  theorem  below,  we  describe  a polyhedral  set  with  a least 
element  and  show  how  this  element  can  easily  generate  a solution  to  the 
problem  (q,  M)  in  the  case  ^ 


- J 
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Theoiem  . I I . IX- fine 


? = { u c : u < 0,  Xu  ^ 0 , Vu  > 0,  c_, u , > - 1 } . 

i=  1 


I 1 


Tfien  7 has  a least  element  u.  Furthermore,  if  then  c u < -1 

u 


and  the  vector  x = -Xi 


T- 
c u 


solves  the  linear  complementarity  problem  (q,  M). 


Proof:  Proposition  2.9  implies  that  ? is  a meet  seini-sublattice  of  R . 
It  is  obviously  closed  and  nonempty  because  o t O'.  We  show  that  O' 
is  bounded  below.  Let  u c clearly  u.  > -l/c.  for  every  i < I 
where  I = {i  : c.  > O}.  Furthermore,  letting  J = fi  : c,  = O},  wo  have, 


(2.  4a) 


ana 


12.  4b) 


’^Jjb  - 'Vi  - ■’^Ji‘^1 


where  d = (d.)  is  the  vector  with  d.  = -l/c,  for  i « I and  o other- 
1 11 

wise.  On  the  other  hand,  it  follows  from  condition  (M3)  that 

T T T T 

i.e.  the  Z-matrices  X^^  and  satisfy  condition  (C2) . Therefore  by 


Lemma  1.2,  the  matrix 


,TJ 

'll 


contains  a complementary  submatrix  which 


is  Minkowski.  This  fact,  together  with  inequalities  (2.4a)  and  (2.4b) 
implies  that  u^  is  bounded  below  by  some  constant  vector.  Hence  o is 
bounded  below.  Therefore  it  has  a least  element  u by  Proposition  2.  lO. 
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Clearly  -u  t 3 If  u f *L( . Thus,  u < -u  for  ever  y u c . Hence 

T—  u 

c u < -1  If  y (}).  It  Is  then  cleo'  that  < V • Therefore  the  last 

c u 

assertaion  of  the  theorem  is  an  immediate  consequence  of  Propositiorr  2.5.  ij 

Remark.  The  assumption  ^ 4>  is  essential  in  order  for  the  vector  x 

T- 

in  the  theorem  to  be  well-defined,  because  otherv/ise,  c u might  be 
zero  (e.g.  c = 0). 

Example  2.12.  Consider  the  problem  (q,  M)  in  Example  2.2.  We  have 
V.  = {(°)}  and 

u,  2 

3'={(^)«R  :u  = 0 and  -2  < u < o}  . 
u.,  1 - 2 - 


/ 0, 

The  least  element  of  J is  (2'* 


This  is  the  solution 
obtained  in  Theorem  2.11 
It  is  also  the  unique  solu 
tion  to  the  linear  program 

(p,q,  M)- 
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which  is  a contradiction.  □ 

In  fact,  the  converse  of  Lemma  2.  13  holds,  namely,  if  (2.  5)  holds, 
then  Ia  = Now,  by  Proposition  1.1,  condition  (2.5)  is  an  equivalent 

formulation  of  condition  (C2).  Therefore,  adding  the  assumption  ^(  = 4> 
to  (Ml)  - (M3)  is  actually  equivalent  to  replacing  conditions  (M2)  and  (M3) 
by  the  single  stronger  condition  (C2). 
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Coroll^iry  2.11.  Let  = «;>.  Then 


T 


u / 0,  Xu  = 0 =>  c u / 0 . 

Proof:  Indeed,  if  there  exists  u nonzero  vector  u such  that  Xu  = 0 
T 

and  c u = 0.  Then  it  follows  from  (Ml)  that  Yu  - 0.  But  this  is 

X 

impossible  because  by  Lemma  2.  13  above  and  Lemma  1.  2,  the  matrix  (^) 

contains  a nonsingular  complementary  submatrix;  this  latter  fact  implies 

u = 0,  contradicting  the  assumption.  □ 

We  have  given  an  example  earlier  illustrating  that  there  can  be 

singular  matrix  X ' satisfying  conditions  (Ml)  - (M3).  The  proposition 

below  shows  that  if  q),  such  matrices  have  rank  at  least  n - 1. 

Proposition  2.15.  If  "U  = q>,  then  rank  (X)  > n - 1. 

Proof:  Suppose  rank  (X)  <n  - 1.  Let  u^  and  u^  be  two  linearly 

1 2 T 1 T 2 

independent  vectors  such  that  Xu  --Xu  =0  and  c u = c u - 1. 

The  existence  of  these  vectors  follows  from  the  assumption  and  Corollary  2.  14. 
We  then  have 

1 T 1 ? T 2 

Yu  + q(c  u ) = Yu  + q(c  u ) = 0 . 


Thus, 


But  we  also  have 


Y(u^  - u^)  = 0 . 


X(u^  - u^)  = 0 


12  12 
Therefore,  u = u . This  contradicts  the  linear  independence  of  u and  u . 

It  is  clear  that  there  are  instances  when  X is  nonsingular,  e.g. 

c - 0.  The  example  below  indicates  that  the  case  where  rank  (X)  - n - 1 
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ts  also  possible  even  If  V • bater  In  our  discussion,  v/e  shall  provide 
necessary  and  sufficient  conditions  for  X to  be  nonsinyular  (sec  Theorem  2.22). 
Example  2.  1 6.  Let  q = (_°)  and  M = (_^  J).  The  problem  (q,  M)  has 

a solution  (°).  Nov/  let  X = (J  ”),  Y = (°  ^),  c = e,  r = (J)  and 

0 T T T 

s = (j^).  Then  MX  = Y + qc  and  r X + s Y > 0.  Therefore  V = 

Nevertheless,  X is  of  ranl<  1.  □ 

In  the  sequel,  we  shall  first  establish  the  lea.st-element  results 
for  the  case  where  X is  nonsingular.  Then  we  shall  prove  two  characteriza- 
tion theorems  having  to  do  with  the  nonsingularity  of  X.  Finally,  we 
shall  investigate  the  last  case,  namely,  where  ‘X  has  rank  n - 1. 

If  X is  nonsingular,  then  the  linear  complementarity  problem  (q,  M) 
is  equivalent  to  the  following  complementarity  problem  (LECP^): 

q(  1 + c u)  + Yu  > 0 , Xu  > 0 and  [ q(  1 + c'^u)  + Yu  ] ^Xu  = 0 . 

The  equivalence  is  based  on  the  transformation  x = Xu.  It  is  clear  that 
X is  a (feasible)  solution  to  the  problem  (q,  M)  if  and  only  if  u = X~'x 
is  a (feasible)  solution  to  the  problem  (LECP^).  Let  3^  denote  the 
feasible  set  of  the  latter  problem,  i.e. 

3^  = {u  € : q(l  + c^u)  + Yu  > 0,  Xu  > O}  . 

T 

Proposition  2.17.  If  V = 4>,  then  1 + c u 0 for  every  u t 3^. 

T 

Proof:  Suppose  there  is  a u < 3^  such  that  1 + c u = 0.  Then  it  follows 

X 

that  Yu  > 0,  Xu^O,  Since  the  matrix  (y)  contains  a Minkowski  comple- 
mentary submatrix,  we  must  have  u > 0 contradicting  the  assumption  that 


T 

1 + c u = 0. 


□ 


By  the  convexity  of  3?^  It  follows  that  one  and  only  one  ot  the 
following  two  statements  must  hold: 


(2.  6a) 
(2.6b) 


1 + c u > 0 for  every  u « 3 


1 + c u < 0 for  every  u « 3^  . 


Lemma  2.  18-  Let  V = 4>  and  3^  ^ <)>•  Suppose  (2.  6a)  ((2. 6b))  holds.  Let 

= {v  € r"  : q + Yv  > (<)  0,  Xv  > (<)  O)  . 

— T 

Then  has  a least  (greatest)  element  v satisfying  c v (>)  i. 

T 

Furthermore,  (q  + Yv)  (Xv)  = 0. 

Proof:  Proposition  2.9  implies  that  is  a meet  (join)  semi-sublattice 

of  R*^.  It  is  obviously  closed.  The  fact  that  (y)  contains  a Minkov/ski 

complementary  submatrix  implies  that  is  bounded  below  (above).  We 

T . , 

show  that  it  is  nonempty.  Let  u « 3^  satisfy  1 + c u > f<)0,  then 
it  is  not  hard  to  see  that  the  vector 


(2.7) 


V = 


1 + c u 


belongs  to  . By  Proposition  2.  10,  has  a least  (greatest)  element  v. 

That  V satisfies  c v < (>)  1 is  clear  because  v < (^  v and  c v < (>)  1 
where  v is  the  vector  defined  in  (2.7).  To  show  that  v satisfies  the 
complementarity  property,  we  refer  the  reader  to  the  proof  of  Lemma  3.  10 
in  [2].  □ 

Remark.  The  polyhedral  set  does  not  depend  on  the  vector  c.  If 

(2.6a)  holds,  this  set  is  precisely  the  set  V mentioned  in  the  introduction. 
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The  followiny  result  cjcncralizos  Theorem  1 . 1 pertuiniiiy  to  tlie 


problem  (q,  M)  with  M t C which  corresponds  to  the  cuse  c - 0. 
Theorem  2,  19.  Let  ~ ({).  Suppose  X is  nonsingular  and  the  problem 

(q,  M)  is  feasible.  Then  the  vector 

(2.8)  X = xl 


VI  - c vl 

where  v is  the  vector  obtained  in  Lemma  2.  l8,  is  a solution  to  tfie 
linear  complementarity  problem  (q,M). 

Proof:  In  fact,  if  x c X(q,  M),  then  u = X c 3^.  Therefore  the 

V 

vector  V is  well-defined  and  so  is  u = ip*  . It  is  easy  to  show 

1 - c V 

that  u solves  the  problem  (LECP^).  Therefore  by  the  equivalence  of 
the  problems  (q,M)  and  (LECP^)  mentioned  earlier,  it  follows  that  the 
vector  X defined  in  (2.8)  solves  the  problem  (q,M).  □ 

Theorem  2.19  shows  how,  under  the  assumptions  in  the  theorem,  a 
solution  to  the  linear  complementarity  problem  (q,  M)  can  be  generated 
from  the  least  (or  greatest)  element  of  the  polyhedral  set  J>^.  One 
notices  that  in  some  instances,  the  linear  complementarity  problem  has  a 
solution  generated  from  the  greatest  element  of  a polyhedral  set.  This 
m.ay  seem,  at  first  glance,  so.mewhat  inconsistent  with  the  title  "A  least- 
element  theory  ..."  of  the  report  and  with  the  phrases  "least-element 
solutions",  "least-element  characterization ",  etc.  ...  which  have  been 
used  throughout  the  report.  However,  it  is  not  hard  to  see  that  such 
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grealost-eleinont  results  can  always  bo  chanyec'  Into  least  -clcmont  results 

by  a very  trivial  modification,  namely,  by  considering  the  negative  of 

the  polyhedral  sets  having  the  greatest  elements.  We  have  not  made  this 

change  of  variables  in  Lemma  2.18,  and  will  not  do  so  later  because  we 

want  to  present  the  results  in  their  most  natural  format. 

In  the  sequel,  we  establish  a relationship  between  the  vector  x 

generated  in  Theorem  2.  19  and  the  solution(s)  to  the  linear  program 

T - 

(p,  q,M)  with  p = r + M s.  Specifically,  we  show  that  x solves  any 
such  linear  program;  moreover,  it  is  the  unique  solution  if  the  vectors  r 
and  s satisfy  the  stronger  condition  (C2). 

If  X is  nonsingular,  then  under  the  nonsingular  transformation, 

X = Xu,  the  linear  program  (p,q,  M)  is  equivalent  to 

T 

(2.9)  minimize  (p  X)u  subject  to  q + MXu  > 0 and  Xu  _>  0 . 


Noting  that 


rri  fTi 

p X = (r  + s M)X  = (r  X + s Y)  + (s  q)c 


we  may  write  (2.9)  as 

T T T T 

(2,  10)  minimize  (r  X + s Y)u  + s q(l  + c u)  subject  to  u € 3^ 


In  order  to  simplify  the  following  discussion,  we  assume  that  (2. 6a) 
holds.  Since  v is  the  least  element  of  it  solves  the  linear  program 

T 

(2.11)  minimize  f v subject  to  v « 


for  any  nonnegative  vector  f,  and  it  is  the  unique  solution  to  (2.  1 l)  for 
any  positive  vector  f.  It  is  then  easy  to  deduce  that  the  vector 


(2.  12) 

solves  the  (fractional)  proyrain: 
(2.13)  minimize 


u = (or  equivalently,  v = ~ ) 

1 - c V 


.3L 
1 -I  c ' u 


f u 


~ subject  to  u < 


1 + c u 

for  any  nonnegative  vector  f,  and  it  is  the  unique  solution  to  (2.  1 3)  for 
any  positive  vector  f. 

Theorem  2.20.  Let  'iq  = 4>.  Suppose  X is  nonsingulai  and  (q.  M)  is 

feasible.  Then  the  vector  x (defined  in  (2.8))  solves  the  lin^'or  program 

T 

(p,  q,M)  with  p - r + M s.  Moreover,  it  is  the  unique  solution  to  any 

T T 

such  linear  program  with  r X + s Y > 0. 

Remark . This  theorem  is  valid  no  matter  v/hich  one  of  t.he  tv/o  inequalities 
(2.6a)  and  (2.6b)  is  valid.  To  simplify  the  proof  below,  we  continue  to 
assume  (2.  6a)  holds. 

Proof  of  the  theorem:  It  suffices  to  show  that  the  vector  u (defined  in 

(2.  12))  solves  (2.10)  and  u is  the  unique  solution  to  (2.  10)  if 

r^X  + s^Y  > 0.  Let  u e 5 , then  t • Thus 

^ 1 + c u 


(r^X  + s^Y)  — - 


T- 
1 + c u 


T T 
< (r  X . s Y) 


T 

1 + c u 


which  implies 


T- 

T T_  T T T T.  T,  T..1+CU 

(2.14)  (r^X  + s Y)u  + s q(l  + c u)£[(r  Xfs  Y)u+s  q(l  + c u)]  ^ 

1 + c u 


Furthormore,  wc  have 


Hence, 


- 11 . jj 

V - . 

1 + c u 1 f c u 


T_ 

< 1 . 

1 + c^u  “ 


On  the  other  hand, 


(r'^X  + s"^Y)  u + s'^q(l  + c'^u)  = r'^(Xu)  + s"^(q  (Y  + qc‘)u>  > 0 . 


Therefore,  for  u « 5^, 

n’T  T T TT  T T 

(r  X + s Y)u  + s q(l  + c u)  < (r  X + s Y)u  + s q(l  -l  c u) 

i.e.  u solves  (2. lO). 

Conversely,  if  u solves  (2.10),  then  (2.  14)  must  hold  as  equality. 
Therefore,  we  have 


(r’^X  + s^Y) ^ = (r’^X  + s^Y) 

1 + c u 


u 

T 

1 + c u 


T T 

i.e.  u also  solves  (2.  13).  Now  if  r X + s Y > 0,  it  follows  from  the 

uniqueness  of  u that  u = u.  [3 

Example  2.21.  Let  q = ( |)  and  M = ( ^ J).  The  problem  (q,  M)  has 

solutions  (J)  and  (°).  If  X = ( ^ 2^’  Lj  4)-  = (J) 

0 T' 

and  c = (^),  then  MX  = Y + qc  and  rX+sY>0.  It  is  not  hard  to 

show  that  M I C . 
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Before  proceeding  to  study  the  ease  lank(X)  - n - ! , we  prove 
two  theorems.  The  first  of  which  characterizes  the  nonsingularity  of  X; 
whereas  the  second  provides  necessary  and  sufficient  conditions  for 
(2.  6a)  and  (2.  6b)  to  hold.  As  a consequence  of  the  theorems,  we  deduce 
that  the  assumptions  in  Lemma  2.  l8  and  Theorem  2.  19  are  actually 
equivalent. 

Theorem  2.22.  The  following  two  conditions  (2.  1 5)  and  (2.  1 6)  ar<^  equivalent: 
(2.16)  X is  nonsingular 

(2.  1 6)  there  exist  a principal  rearrangement,  with  permutation  .matrix  1, 
and  partitioning  of  X,  Y,  q and  c such  that 


(2.  I6a) 


T 

P XP  = 


11 

T 

Y 

12  ' 

X i 

, P YP  = 1 

lY 

Y.,.,  i 

21 

221 

\ 21 

221 

(2.  I6b) 
(2.1 6c) 


Xii  and  (^22  - are  nonsingular 


(2.  16cl) 

Proof:  That  (2.  16)  implies  (2.  l6)  is  trivial.  For  the  converse,  it  suffices 
(see  [1])  to  show  that  the  matrix  ^22  ' ^21^11^12  >^°'^singular.  We  have 


“l2\ 

'"n 

"12I 

^.2\ 

“21 

“22/ 

\’'21 

i'^2. 

’'22. 

T 

T I P 
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Solving  for  In 


“2/.,  * “2/2,  - ’'2,  ^ V. 


gives 


= (Y^  - M,,X,  )X"^  +q,c’^X"‘  . 

21  ' 21  22  21'  11  ^2  1 11 


Substituting  this  latter  equality  in 


"2/12  ^ >^22^22  = ’'22  * ‘>2'^2 


yields 


“22<>‘22  - = <’^22  - ’^21>';'.Xl2'  ^ ‘'2<‘=I  ‘ =X^.2' 

' <’'22  - ’'2.^;^2'l  ‘ " •''22  - ’'2l’‘n\2'“‘’2'=I  - =X!^2>1  • 

Condition  (2.  l6d)  implies 

detll  MV,2  - V,^x;;x^2>' VI-  =XK2>1 
= >*<"I-=IVl2><’'22-  Vu=‘l2>''2'°- 

This  establishes  the  nonsingularity  of  X^^  - ^21^11^12’  ^ 

Remark.  The  argument  used  above  is  a generalization  of  the  one  in 
Theorem  1.3  for  the  case  M c C . In  the  proof,  we  have  used  the 
condition  (Ml)  only.  So  in  fact,  the  theorem  provides  necessary  and 
sufficient  conditions  for  the  matrix  X (not  necessarily  Z-matrix)  satisfying 
MX  = Y + qc  for  some  matrices  M and  Y (not  necessarily  Z-matrix) 
and  some  vectors  q and  c (not  necessarily  non-negative),  to  be 
nonsingular. 
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L 

r 

[ Theorem  2.23.  Let  y - 4'.  Suppose  3^  ^ 4»-  Then  {l.hu)  holds 

i if  and  only  If  there  exist  a principal  rearrangement  with  permutallon  matrix 

P,  and  partitioning  of  X,  Y,  q and  c such  that 


X 

X \ 

fY 

Y \ 

(2. 17a)  P^XP  = 

I 1 
X 

12) 

X 

T 

, P YP  r 

11  12 
Y,  Y, , 

^ 21 

zzi 

V 21  22/ 

T 

1 nT 

(2.17b) 

P q = 

g 

1 , P c = 

c 

1^2 

1 2 

(2.17c) 


(2.  I7d) 


> (^)O 


Proof;  "Sufficiency".  Let  u < 3^.  Then  we  have 


q^d  + c u)j 


Condition  (2.17c)  implies  that 


Therefore  1 + c^u  >(<)0  for  every  u « 3^,  i.e.  (2.6a)  ((2.6b))  holds. 
"Necessity”.  The  existence  of  the  permutation  and  partitioning  such  that 
conditions  (2.17a)  - (2.  17c)  are  satisfied  is  an  immediate  consequence 
of  the  assumption  “U  = <t>  and  Lemma  1.2.  So  it  remains  to  verify  (2.  I7d). 
The  deduction  above  shows  that  (2.20)  is  valid  for  every  u * 3^.  Now 
if  3j  ^ <|>  and  (2.6a)  ((2.6b))  holds,  then  it  follows  from  (2.20)  that 
1 + c^t  > (<)0,  or  equivalently,  (2. 17d)  holds.  □ 

We  have  seen  that  if  X is  nonsingular  and  if  the  problem  (q,  M) 
is  feasible,  then  3^^  * <t>.  Combining  Theorems  2.22  and  2.23,  we 
conclude  that  if  V = <l>  and  if  3^  ^ 4>  then  (q,  M)  is  feasible  (in  fact, 
has  a solution)  and  X must  be  nonsingular. 


_ — n 

Corollary  2.24.  If  “K  - 4>,  the  following  arc  equivalent 

(2.21)  X Is  nonsingular  and  (q,  M)  has  a solution  (given  by  (2.  8)) 

(2.22)  3^  <|>  • 

2.4.  The  case  where  V = <i>  and  X is  singular.  In  the  rest  of  this 

paper,  we  investigate  the  remaining  case,  namely  = 4>  and  rank  (X) 

= n - 1 . Let  u^  be  the  (unique)  vector  satisfying 

0 TO 

Xu  =0  and  c u = 1 . 

The  existence  of  such  a vector  follows  from  Corollary  2.  i4.  Then,  it 
follows  that 

(2.  23)  q = -Yu°  . 

1 T 1 

Let  u be  a nonzero  vector  satisfying  Xu  =0.  It  follows  from 

elementary  linear  algebra  that,  for  any  x € r’’,  there  exist  u < 

and  « € R^  such  that 

(2.  24)  X = Xu  + au^  . 

Under  such  a representation,  we  have,  by  (2.  23) 

q + Mx  = q + MXu  + aMu^ 

T 1 

= q(l  + c u)  + Yu  + aMu 

TO  1 

= Y[  u - (1  + c u)u  ] + aMu  . 

Proposition  2.25.  Under  the  identification  (2. 24),  the  linear  complementarity 
problem  (q,  M)  is  equivalent  to  the  problem  (LECP^)  of  finding  a vector 
u t and  a scalar  a such  that 
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aMu^  i Y[  u - (1  + c*^u)u^]  > 0 
+ X[  u - (1  + c '^u)u*^  J >0 

[ oMu^  + Y(u  - (1  + c ^u)u^)  ] ■*[  ■(  X(u  - (]  f c "u)u^)  ] = 0 . 

Let  3^  denote  the  feasible  set  ot  (LLCP^),  i.e. 

5-  = {(^)  < : tvMu^  + Y(u  - (1  + c^u)u^)  > 0,  au^  +X(u  - (1  + c^u)u^)  > 0 } 

Z a 

Proposition  2.26.  If  h = <J),  then  a * 0 for  every  (^)  t 32- 

Proof:  Suppose  there  exists  (^)  « 3.,  such  that  a = 0.  Then 

a Z 

Y[  u - (1  + c^u)u^  ] > 0 
X[  u - (1  + c^u)u*^]  > 0 . 

T 0 

Hence,  by  Lemma  1.2,  it  follows  that  u - (1  + c u)u  > 0.  Thus 

T TO 

0<c  [u-(l+c  u)u  ] = -1 

which  is  a contradiction.  □ 

Remark.  An  equivalent  formulation  of  the  conclusion  in  Proposition  2.26 
can  be  stated  as  X(q,  M)  D R(X)  = c}). 

By  the  convexity  of  it  follows  that  one  and  only  one  of  the 

following  two  statements  must  hold: 

(2.  2 5a)  o > 0 for  every  (^)  « 3^ 

(2.  25b)  a < 0 for  every  (^)  « 3,  . 

a c 

Lemma  2.  27.  Let  = 4>  and  Suppose  ( 2.  25a)  ((2 . 25b))  holds 

Define 

{v  t r'^  : Mu^  + Yv  > (<)  0,  u^  + Xv  > (<)  0}  . 
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Then  has  a least  (greatest)  element  v satisfying  c v<(>)0. 

Furthermore,  (Mu^  + Vv)^(u^  t Xv)  = 0. 

1 T 0 T 1 

Proof;  If  (^)  € then  v = ~(u  - (1  + c u)u  ) t and  c v = - 

The  rest  of  proof  proceeds  in  the  same  way  as  in  the  proof  of  Lemma  2.18. 

Theorem  2. 28  . Let  = <}>.  Suppose  rank(X)  = n - 1 and  the  problem 

(q,  M)  is  feasible.  Then  the  vector 

(2.26)  X = - ~^(>^  + u^) 

c V 

where  v is  the  vector  obtained  in  Lemma  2.27,  is  a solution  to  the  linear 
complementarity  problem  (q,M). 

Proof:  The  assumptions  in  Lemma  2.27  are  satisfied.  Therefore,  x is 

V ^ 1 

well-defined.  It  suffices  to  show  that  u = - and  a = - solves 

c V c V 

the  problem  (LECP^).  Noting  that  1 + c u = 0,  we  can  easily  deduce 
u 

that  (_)  indeed  solves  (LECP_).  This  establishes  the  theorem*.  □ 

a 2 

Having  demonstrated  how,  under  the  assumptions  of  Theorem  2.28, 

a solution  to  the  linear  complementarity  problem  can  be  generated  from  the 

least  (or  greatest)  element  of  the  polyhedral  set  we  next  establish 

a relationship  between  this  solution  x and  the  solution(s)  to  the 

T 

linear  program  (p,  q,  M)  with  p = r + M s. 

Theorem  2.29.  Let  V = <!>•  Suppose  rank  (X)  = n - 1 and  (q , M)  is 

feasible.  Then  the  vector  x (defined  in  (2.26))  solves  the  linear  program 

T 

(p,  q,M)  with  p = r + M s.  Moreover  it  is  the  unique  solution  to  any 

T T 

such  linear  program  with  r X + s Y > 0. 
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r.- 


Proof;  Assume  for  simplicity,  that  (2.25a)  holcl.s.  The  aryument  which 

follows  is  similar  to  that  of  Theorem  2.  20.  Let  x = Xu  + au  t X(q,  M), 
1 T 0 

then  V = ~(u  - (1  + c u)u  ) i 3 ,.  Hence  we  have 
a c 

_ 1 TO 

V < “(u  - (1  f c u)u  ) 

~ a 

which  implies 

T_  1 
c V < - — . 

~ a 


Moreover,  by  (2.23), 

(2.27)  (r'^X  + s'^Y)v  < ^(r^X  + s'^Y)(u  - (1  + c'^u)u*^j 

ITT  T T 

= -[(r  X + s Y)u  + s q(l  + c u)j  . 

ot 

Furthermore, 

r'^(Xu  f au^)  + s^[q  + M(Xu  + au*^)]  > 0 . 

Combining  the  above  inequalities,  we  deduce, 

p^x  = (r^  + s"^M)(-  -:^(X7  + u^) 

C V 

= - + s^Y)v  f s'^q(c'^v)  + (r'^  + s'^M)u^l 

C V 

< - ~~[;~{(r'^X  + s^Y)u  + o'(r'^  + s'^M)u^  + s'^q(l  + c u)}  ] - s q 

T—  ^ 
c V 

= (-  -^)(-)[r^(Xu  + aub  + s’^(q  + M(Xu  + au*))]  - s^q 

C V 

T T T 

< (r  + s M)x  = p X . 

Therefore  x solves  the  linear  program  (p,q,M). 
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Conversely,  if  x Is  a solution  to  (p,  q,M),  then  we  must  have 

T T - 

equality  in  (2.27).  If  r X + s Y > 0,  then  by  tfio  uniquonf-ss  of  v,  it 


1 


T . 0, 


follows  that  V = “(u  - (1  + c u)u  ).  Thus  a - - -tt'  and  Xu  = - Xv. 


Therefore,  x = Xu  + au  = - “~(Xv  + u ) = x.  This  establishes  the 

c V 

theorem . o 

We  conclude  this  paper  by  proving  a result  (Corollary  2.  31)  which 
strengthens  Corollary  T of  Mangasarian's  report  [?]  where  it  was  shov/n 
that  if  n > 3 and  M.  is  a positive  matrix  which  is  diagonally  dominant 
colutnn-by-column,  that  is,  m..  > Y.  m,.,  j = 1,  . . . ,n,  then  conditions 

’ ij’ 

(Ml)  - (M3)  are  satisfied  (in  fact  for  every  n-vector  q). 

Proposition  2.  30.  Let  M be  an  n x n matrix  satisfying 


(2. 28a) 


MX  = Y for  some  X,  Y c Z 


(2.  28b) 


r X > 0 for  some  r > 0 


(2.  28c) 


Y = (Y  ) is  such  that  Y. . < 0 for  every  j ^ i 
ij  iJ 


Then  M « C. 

Proof:  Let  e > 0 be  small  enough  such  that  Y + eM  < Z.  Then 
M(X  + el)  = Y + eM  and  X + el  « K.  Therefore  M < C.  □ 

Corollary  2.  31.  Let  n > 3 and  let  M be  a positive  n x n matrix 
which  satisfies  either  of  the  following  conditions: 

(2.29a)  diagonal  dominance  column-by-column:  m,.>  ^m  , j = l,...,n 

JJ  iJ 

(2.29b)  diagonal  dominance  row-by-row:  m > Y ^ " l,...,n  . 


Then  M t C. 
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BB— ~ 


Proof;  Wo  may  without  loss  of  cjcMioraiily,  assume.-  Uiat  M has  been 

normalized  so  that  it  has  ones  on  the  diagonal.  Lot  M I ) F where 

f = (L,  ).  = 0 and  F,.  = m,.  for  i ^ j.  Define  X - I - F and 

Ij  li  1)  1} 

Y = MX  = I - F^.  Clearly  X and  Y « Z.  Furthermore,  it  is  easy  to 
see  that  condition  (2.28c)  is  satisfied.  We  show  that  (2.28b)  is  also 
satisfied.  It  is  clearly  satisfied  if  (2.29a)  holds.  Now,  if  (2.29b)  is  true, 
then  X € Kq  (by  Theorem  (5. 4)  in  [ 3])  and  X is  irreducible  (see  [ 8 ]) . 
These  two  properties  of  X imply  (2.28b)  (by  Theorem  (5.8)  in  [ 3]). 
Consequently,  by  Proposition  2.30,  M c c.  □ 
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